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As a prototype model of antiferromagnetism, we propose a repulsive Hubbard Hamil- 
tonian defined on a graph A = AVJB with AC\B — and bonds connecting any element 
of A with all the elements of B. Since all the hopping matrix elements associated with 
each bond are equal, the model is invariant under an arbitrary permutation of the 
^-sites and/or of the S-sites. This is the Hubbard model defined on the so called 
{Na, Ai's)-complete-bipartite graph, Na (Nb) being the number of elements in A (B). 
In this paper we analytically find the exact ground state for Na = Nb — N a.t half fill- 
ing for any A'^; the ground state expectation value of the repulsion term has a maximum 
at a critical A'^-dependent value of the on-site Hubbard U and then drops like 1/U for 
large U. The wave function and the energy of the unique, singlet ground state assume 
a particularly elegant form for A' — > oo. We also calculate the spin-spin correlation 
function and show that the ground state exhibits an antiferromagnetic order for any 
non-zero U even in the thermodynamic limit. This is the first explicit analytic example 
of an antiferromagnetic ground state in a Hubbard-like model of itinerant electrons. 
The kinetic term induces non-trivial correlations among the particles and an antiparal- 
lel spin configuration in the two sublattices comes to be energetically favoured at zero 
Temperature. On the other hand, if the thermodynamic limit is taken and then zero 
Temperature is approached, a paramagnetic behavior results. The thermodynamic 
limit does not commute with the zero- Temperature limit, and this fact can be made 
explicit by the analytic solutions. 



I. INTRODUCTION 

The Hubbard Hamiltonian is one of the most popular models to describe strongly correlated electron 
systems. In spite of its simple definition, it can be exactly solved only in few cases. An example is the 
milestone solution of the Hubbard ring [Q] by means of the Bethe-^nsate extended to fermions |^ [^. 
Other exact solutions were worked out for the one-dimensional Hubbard ring in the presence of an external 
magnetic field for the one-dimensional Hubbard chain with open boundary conditions Q and thereafter 
for the SU{N) one-dimensional Hubbard ring Q. However, when the space dimensionality becomes bigger 
than 1, few exact results are available and usually they concern the ground state properties. Among them, 
we mention the Lieb theorem on the ground state spin degeneracy fel that will be explicitly used in this 
work. Exact ground state wave functions are even more infrequent. To the best of our knowledge the only 
non-trivial results are the ferromagnetic ground-state solutions devised by Nagaoka (infinite repulsion and 
one hole over the half filled system), by Mielke and Tasaki jlo| [QtI (for which the kinetic energy spectrum is 
macroscopically degenerate or at least very flat) and by Wang (infinite-range hopping). 

In the Hubbard Hamiltonian with infinite-range hopping a particle in any site can hop to any other site of 
the system; the associated graph is said to be complete. This model was numerically studied by Patterson on 
small clusters in 1972 and solved in the thermodynamic limit by van Dongen and VoUhardt only at the 
end of the eighties. Much more effort was needed to find the exact ground state(s) in the finite-size system. 
Verges et al. managed to accomplish this task for arbitrary numbers of particles and sites in the limit of 
infinite on-site repulsion U by exploiting a scheme proposed by Brandt and Giesekus Two years later 
Wang constructed explicitly the ground states of the system for arbitrary U and number of particles above 
half filling In the case of one particle added over the half-filled system the ferromagnetic ground-state 

solution follows as a special (and the easiest) case of general results by Mielke and Tasaki [ p^ . 

In this paper we find the exact ground state wave function of the half-filled Hubbard model on the Complete- 
Bipartite Graph (CBG) for arbitrary but finite U. The CBG K = AVJB has bonds connecting any element of 
A with all the elements of B and can be considered as the natural further step (with respect to the complete 
graph previously described) towards the standard Hubbard Hamiltonian (defined on the hypercubic lattice and 
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hopping between nearest-neighbours sites). Even if the CBG is still somewhat unrealistic, we feel that exact 
solutions should always be welcome, expecially because they lend themselves to be generalized. Furthermore, 
our solution is the first example of antiferromagnetic ground state in a model of itinerant electrons, in contrast 
with the ferromagnetic solutions mentioned above; it may provide useful hints about antiferromagnetism 
outside of the strong coupling regime (where the Hubbard model can be mapped onto the Heisenberg model). 

The paper has been organized as follows. In Section |^ we define the model and discuss the physics of 
the non-interacting (U = 0) Hamiltonian together with few relevant examples of finite-size realizations. In 
Section III we study the thermodynamic limit of a class of Hubbard-like models, which includes the Hubbard 
Hamiltonian on the complete graph and on the complete bipartite graph, having a non-extensive number of 
isolated one-particle energy levels plus a single level whose degeneracy is proportional to the size of the system. 
Following the reasoning of van Dongen and VoUhardt , we show that the kinetic term is totally decoupled 
from the interaction term and that the system behaves as a paramagnet: the spin-spin correlation length is 
zero. These results are a consequence of the trivial behaviour of such models any time the thermodynamic limit 
is taken first. On the contrary, much more difficult is to find exact properties in finite-size systems and diff'erent 
answers may be obtained if the thermodynamic limit is taken only at the end of the calculations (as we shall 
show in Section |v^ . In Section ^ we introduce the essential tools to face the problem. Let N = Na ~ Nb 
be the number of sites of each sublattice A and B. First, we find a (2A'^ — 2)-body determinantal eigenstate 
|$>if) of the Hamiltonian with vanishing double occupation; then, we demonstrate that it is a key tool to 
build the ground state. We show that mapping the ^-sites onto the Z3-sites and viceversa, \$af) retains its 
form except for a spin-flip; we shall call this property the antiferromagnetic property. Several analogies with 
the properties of the half-filled standard N x N Hubbard model are pointed out at this stage. Here we also 
deal with the spin projection of |$af) onto the singlet and the triplet subspace and useful identities between 
the two spin-projected states are obtained. Then, in Section ^ we propose an Ansatz for the ground state 
wave function at half filling containing the singlet and triplet projections of \$af)- We set up the Schrodinger 
equation and by exploiting the antiferromagnetic property we close the equations and get 3 exact eigenstates. 
The ground-state uniqueness proved by Lieb Q| is used to show that the lowest-energy state of our Ansatz 
actually corresponds to the ground state of the system. Remarkably, the ground state energy is negative 
for any value of the repulsion U ; qualitatively, we may say that the particles manage to avoid the double 
occupation very effectively. In Section ^ we study the ground state energy Eo as a function of U and of the 
volume of the system and we discuss the implications of exchanging the thermodynamic limit with the limit 
of zero temperture. We find that Eq is a monotonically increasing function of U and A'' due to the existence 
of non-trivial correlations even for large A''. The nature of these correlations is investigated by computing 
the expectation value of the repulsion. We show that for any finite N there is a critical value of U yielding 
maximum repulsion. From the exact spin-spin correlation function we find that the ground state average of 
the staggered magnetization operator squared is exstensive. Hence, the ground state is antiferromagnetically 
ordered; we underline that this holds not only at strong coupling, but for any value of U . Finally, a summary 



of the main results and conclusions are drawn in Section VH 



II. DEFINITION OF THE MODEL 



Let A = ^ U S with ^ n S = be a collection of sites and 

\A\=Na, \B\=Nb ^\h-\^NA + NB ■ (1) 

Here and in the following we shall denote by \S\ the rmmber of elements in the set <S. We consider the Hubbard 
Hamiltonian 

-H'Hubbard = -f^O + W, ^0= '^^tx ,ycl^Cy„ , W = U'^^n^^h^l (2) 

where Cxa {c\cr) is the annihilation (creation) operator of a particle at site x with spin a [ and fixer ~ c^^Cxa 
is the corresponding particle number operator. The hopping matrix T with elements Tx.y — tx.y = ty^x is a 
real-symmetric matrix while U is a. positive constant. If 

^ _ f t ior X & A{x £ B) and y £ B {y £ A) , . 

^'^ I otherwise 

the graph A is said to be complete bipartite and we call the Hamiltonian in Eq.(|2|) the CBG-Hubbard Hamil- 
tonian. The model is invariant under an arbitrary permutation of the ^-sites and/or of the B-sites. Therefore, 
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the symmetry group includes Snj^ (8) Sjv^ , Sn being the set of permutations of A*' objects. As usual, the full 
group is much bigger: the presence of spin and pseudospin symmetries leads to an 5*0(4) internal 

symmetry Group |l9j and in the case Na = Nb there is a Z2 symmetry because of the A ^ B exchange. 

In Fig.(|l]) we have drawn a few examples of finite-size systems. For A'^ = 1 and N = 2 the model is equivalent 
to a one dimensional ring of length L = 2, 4 respectively. For A'' = 3 we have what can be understood as a 
prototype, (1,1) nanotube model, the one of smallest length L = 1, with periodic boundary conditions. For 
general A'', one can conceive a gedanken device, like the one illustrated pictorially for N = 6 in Fig.(|l|.d). 
(The whole device should actually have the topology of a torus, and the two horizontal faces should coincide, 
but this is not shown in the Figure for the sake of clarity.) The A'^ vertical lines represent a realization of 
the A sublattice while the B sublattice is mimicked by the central object. The radial tracks in the Figure 
represent conducting paths linking each A site to each B site according to the topology of our model. The A 
sites are represented by one-dimensional electron-in-a-box systems of length L. Each A site has the one-body 
energy spectrum efroc ~ /mL^ where n is an integer and m is the electron mass. We assume that L is so 
small that the excited states are at high energy and can be disregarded in the low-energy sector; this requires 
mL^ < 1/(7. Here, of course, U denotes the Coulomb self-energy of a box with two electrons. The B sites 
are hosted by A'^ quantum dots, that can be represented by 5-function-like attractive potential wells of depth 
V; if \V\ is large and the radius of the dots is ^ L, these are practically independent of each other. The 
one-body energy of each B site is Ebound ~ —mV^\ the unbound states can be neg lected if mV^ > U. We 
are assuming for simplicity that the U of the B sites is the same as for the A ones. Turning on a constant 
potential Vo on the central object, one can arrange that the Hubbard Hamiltonian on the CBG Hnubbard gives 
a good description of the system, with filling < 2 per site. 



t N=6 




FIG. 1. Physical representation of the CBG for N = 1 (a) and N = 2 (b) which is equivalent to a one dimensional 
ring of length L = 2, 4 respectively. For = 3 (c) we have the (1,1) nanotube of smallest length and periodic boundary 
conditions. For N = 6 (d) we present the gedanken device described in the text. 



It is convenient to label each site with an integer in such a way that x = 1, . . . , Na corresponds to sites in 
A and X = Na -I- 1, . . . , Na + Nb corresponds to sites in B. In the special case Na ~ Nb ~ N the hopping 
matrix can be written as 

where Z is the N x N null matrix and J is the N x N matrix whose generic element J^^y = 1. The one-body 
spectrum has three different eigenvalues 

£g = -Nt = 

£0 = 

Eg = Ni=t (5) 

with degeneracy dg = 1, do — 2N — 2 and dg — 1 respectively. We use the convention t > so that Eg is the 
lowest level and we shall call Shf the set of zero-energy one-body eigenstates. 
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The orthogonal matrix O that diagonahzes T can be written in the form 



O - 



V 



1 1 
R 

_L 1 

/2N V2N 



-1 -1 
/2N \/2N ■ 

R 

_2 1 

/2N \/2N ■ 



(6) 



where Zr is the (TV — 1) x TV rectangular null matrix, while R is an (TV — 1) x TV rectangular matrix whose rows 
are orthonormal vectors which are orthogonal to the TV— dimensional vector (1,1,... , 1, 1). The zero-energy 
one-body orbitals can be visualized by a simple argument. Consider any pair x, y, with x ^ y, of sites 
belonging to the same sublattice, say A, and a wavefunction tf}x,y taking the values 1 and -1 on the pair and 
elsewhere in A and in B. It is evident that tpx,y belongs to Shf- Operating on il)x,y by the permutations of 
Sn we can generate a (non-orthogonal) basis of TV — 1 eigenfunctions[] vanishing in B; further, by means of 
the Z2 symmetry, we obtain the remaining orbitals of Shf, which vanish on A. This exercise shows that the 
group considered above justifies the (2TV — 2)-fold degeneracy of the one-body spectrum. 

Having obtained the one-body basis, one can form a suitable R matrix by orthogonalization. Let us define 
the transformed operators: 



2N 

g = ^ Oi,xCx 



x = l 



E{Cx — Crr-l-jv) 

x^\ 

ai = Oi+i,xCx = Ri,xCx, i — I,. 

a:— 1 x^l 
2N N 

= ^ ^ Oi+N,xCx ~ ^ ^ Ri,xCx+N, i — 1, ■ 



,TV- 1 



x = l 

2N 

g = 02N,x 

x = l 




,TV- 1 



Cx + Cx+n) 



(7) 



The inverse transformation reads 



^=ig + g) + Ri.xcii, xgA 

i 

/2TV " ^ ' 



and the kinetic term Hq becomes 



Ho 



igcr - gtga) 



(8) 



(9) 



Hence, if we do not rescale the hopping constant, the average kinetic energy remains an extensive quantity 
proportional to Ni: the kinetic energy of the two particles in the lowest level coincides with the kinetic energy 
of the whole system. On the other hand, if f ~ 1/TV the two energy gaps eo — £g and Eg — eo remain constant 
as TV increases. As we shall show in the next Section, the model can be exactly solved in the thermodynamic 
limit whatever is the dependence of £g and £g on |A|. 



III. THERMODYNAMIC LIMIT 



The Hubbard Hamiltonian on the CBG (as on the complete graph) belongs to a class of Hubbard-like 
models that can be exactly solved in the thermodynamic limit. Here we generalize the results by van Dongen 



^ There is a simple trick to see that the number of independent wavefunctions obtained in this way is TV — 1. 
Let P £ Sn be a permutation such that P{x) = x and P{y) 7^ x. We get the full basis since P{y) takes TV — 1 
values with independent tpx,p{y)- 
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and VoUhardt to the case of a non-extensive number of non-vanishing one-body eigenenergies and in the 
presence of a local external magnetic field coupled to the spin of the particles. Even if the reasoning is similar 
to that of Ref a detailed presentation of the main results is needed to clarify more subtle points and 
make our work self-contained. Furthermore, some of the results are absent in the original paper and we believe 
that their derivation will facilitate the reader in the comprehension of what follows. 
Let us consider the Hamiltonian 

H — -ffnubbard — M JT-ict " h^in^t — W-xi) = J^Hubbard + + Hh (10) 

where -ffnubbard = ^^o + W is defined in Eq.(j^, is the chemical potential and Hh represents the coupling 
with an external magnetic field. Let us assume that only a finite number f of eigenvalues of the hopping matrix 
T are non-vanishing, £i 7^ for i = 1, . . . ,f, while the remaining |A|— f are identically zero. Denoting with <j>i 
and (f)\ the annihilation and creation operators that diagonalize we have 

f 

^^0=^^6.010.. (11) 

i — 1 CT 

and the CBG corresponds to the case f=2, see Eq.(^). The Gran-Partition Function -Z^[T, U, {h^}] can always 
be written as 

Zp[T,U,{h^}] = Zp[T,0,{h^ = Q}^^y^,[T,u,{n.}]^ Wp[T,U,{h^}] = (e-/o'"^'"''")S (12) 

where 13 is the inverse Temperature, Hiij) = Wij) + Hhir) is the interaction Hamiltonian in the imaginary- 
time Heisenberg representation and (. . .)o is the thermal average in terms of 2/3 [T, 0,{ha,= 0}] = [T] where, 
by the linked-cluster theorem, one has to retain only those contributions represented by connected diagrams. 
As far as f remains fixed, we can safely substitute Vy^[r, U, {/i^}] with W^[r = 0, ?7, {h^Y\ in Eq.(^. Indeed, 

W/5[T,i7,{ft41 = W^[r = 0,(7,{/i4] + O(f/|Al) (13) 

as long as the unperturbed one-body Green's functions Gicr(j^ ~ {T(pia{T)cf>l^{0))o, i = 1, . . . ,f, do not diverge 
in the thermodynamic limit; here T is the Wick time-ordering operator and {. . .)o is the thermal average with 
U — and {hx =0}. From the explicit expression of daij), 

_e-fe-M).j^ r<0 ' ^^^> 

with /i = [1 4- e''^^'"^']"^ the Fermi distribution function, it follows that Giair) converges to a finite value 
whenever lim|A|^ooei is well defined. One can use the result in Eq. (p^) to write the exact Gran- Partition 
Function in the thermodynamic limit 

Zp[T,U,{K}] = ^^^^Z^[T^O,U,{hx}] = 

-ri[ '^%';:/""1 n [l + 2.cosh,..+.^e-n , (15) 



where z — e''^ is the fugacity, and the thermodynamic potential f2^[T, U, {hx}] = —f3 ^ InZ^lT, U, {hx}] 
Qp[T, U, {hx}] = i 1 2f ln(l + z)-Y^ ln[l + 2ze-^^^ + z^e"""-] -^\n[l + 2z cosh/?ft, + z'e"''^] 1 . 

^ i ^ = l xeA J (16) 

Eq.(|l6|) reduces to the result of van Dongen and VoUhardt in the case i=l, hx — £ A and ei oc — |A| < 0. 
All the thermodynamic quantities can be derived from Eq.(^^ and we defer the interested reader to the 
original paper by van Dongen and VoUhardt for further details. Here we calculate the magnetization at site 
X, nix ~ —^{dflfi/dhx) = ^{fixT; — fixi) = (Sx) and the connected spin-spin correlation function {SxSyY = 
(SxSy) — {Sx){Sy) = —{l/4:f3){d^Qi3/dhxdhy), where (. . .) means the average in the gran-canonical ensemble. 
Due to the exact decoupling of the kinetic energy, we expect from the begining a paramagnetic behaviour. 
We have 



1 do, IS z sinh /3hx 



2 dhx 1 + 2z cosh I3hx + z^e-f^^ 



(17) 
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where z can be expressed in terms of the number density by means of the relation 



_ 1 dQ.p _1z^ cosh f3h^ + ze-f^" 
" ^ ~|A| "a^ ^ iAl ^ I + 2z cosh I3h^ + z^e-i^u ' ( > 

Since z is finite for any non-zero Temperature, nix when hx 0, i.e. 

lim lim hm nix = . (19) 

/3— »oo h^~*0 |A|-^oo 

On the other hand, nix 7^ if we exchange the last two limits in Eq.(^9|). Indeed, taking hx = /i=const> 
for the sake of clarity, Eq.(p^ can be exactly solved for z: 

V(n - 1)2 cosh^ Ph + n{2 - n)e-0^' - (1 ^ n) cosh f3h f ^Z^,"^^^"^^ ^^'^,2 

"^"^ [ e-'3'^n/(l ~n) n < 1 (20) 

Substitution of these asymptotic behaviours in Eq.(|l7|) yields 




lim lim lim nrix = < 1/2 n = 1 . (21) 

/ij;^0+ /3->oo |A|-.oo 

The case /i < is similar and one can show that Eq.(^l|) remains true in the limit hx — > 0~ if m,x --nix. 
Therefore, the magnetization does not depend on U and the zero- Temperature limit of the model is the same 
of the paramagnetic Hamiltonian ffpara = + Hu, where and Hh are defined in Eq.(0). From Eq.(0) 
we also conclude that 



a result which is independent of the local external field configuration: in the thermodinamic limit two localized 
spins do not interact and the spin-spin correlation length is zero. Denoting with Gapin(a;, y; /3, h, |A|) the thermal 
average of S^Sy at inverse Temperature /3, external field h and size of the system |A|, from Eq.(|l9|) we get 

lim lim lim (S^Sy) ^ lim lim lim GapiJx,y; (3, h,\A\) ^ , (23) 

that is, the spin-spin correlation function vanishes if we first take the limit h ^ and then the limit of /3 ^ cxj. 
We emphasize that the above results are correct only if we first take the limit |A| 00. The thermodynamic 
limit makes the model trivial and different graph structures, like the complete graph and the complete bipartite 
graph, all yield the same thermodynamic behaviour. On the contrary, a much more difficult task consists in 
finding exact results in finite-size systems. To face such problems is not only a mathematical exercise. After 
that in Section ^ the exact ground state of the half-filled Hubbard model on the CBG is explicitly written 
down for arbitrary value of U and |A| and in the absence of an external field h, we shall see in Section ^ that 
the thermodynamic limit, |A| — > 00, and the limits of zero Temperature and external field, (3 00, h 0, 
do not commute. 



IV. THE W = STATES 



In this Section we introduce the essential tools to deal with the antiferromagnetic ground-state solution. Let 
us consider the one-body eigenstate aj|0) with vanishing amplitudes on the B sublattice, that is (Ojca^aJjO) — 
if X £ B. Similarly, feJjO) has vanishing amplitude on the A sublattice and therefore the {2N — 2)-body state 

|$i^i,) = aL...ajv_i,6L...6ki.|0), ^ = (24) 

is an eigenstate of Ho and of W with vanishing eigenvalue. In the following we shall use the wording W — 
state to denote any eigenstate of //nubbard in the kernel of W . It is worth to observe that by mapping the 
^-sites onto the B-sites and vtceversa, \^af) retains its form except for a spin-flip; we call this property the 
antiferromagnetic property for obvious reasons. 

In the non-interacting {U = 0) half filled system, the structure of the ground state is trivial: two particles 
of opposite spin sit in the lowest energy level g, while \Shf\ ~ (27V — 2) particles lie on the shell Shf of zero 
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( 2N — 2 \ ^ 

energy. In the spin — subspace this ground state is ( ^ ^1 times degenerate. To first order in W , 

the degeneracy is only partially removed^ Indeed, if Ps,Ms is the projection operator onto the subspace of 
total spin S with 2-component Ms, the structure of the determinantal state in Eq.(^) yields 

PsA^^aI) = I-I-af) / 0, V5 = 0, . . . ,iV - 1 ; (25) 

therefore, the states {(7|'?||'&^']!'), S = 0, . . . , — 1} belong to the ground state multiplet in first or- 
der perturbation theory (being = states). Since the Lieb's theorem Q ensures that the interacting 
ground state \'^o{N,U)) is a singlet, only (;|(7j 1$'^'^) can have a non- vanishing overlap with \^q{N,U)) in 
the limit i7 0"*". In Section ^ we shall prove that (?j(/j|"J>^^) is the ground state for U = 0'^ , that is 

lim(/^o{*I'o(-^, f^)l5-[<?||3>^^) — 1. In this model, the first-order solution has a peculiar significance. Usu- 
ally, the exact ground state has no overlap with the non-interacting one, in the thermodynamic limit: this 
is because the (7 = ground states define a proper subspace of the full Hilbert space. Below, we prove that 
{'I'o(A'^, t^)l5|s'||'J'^F) 7^ ' that is, the lowest approximation keeps a finite weight. 

At this stage, we note the analogy of the above results with those relevant to the standard Hubbard model 
defined on a A'^ x A'^ square lattice with periodic boundary conditions (hopping matrix elements only between 
nearest neighbor sites). The determinantal state resembles the W = state that we obtained for even 

N in Ref. |^^. The degeneracy of the zero-energy one-body eigenspace was 2A'^ — 2 and it was shown 
that A*' — 1 zero-energy eigenfunctions have vanishing amplitudes on a sublattice while the remaining A^ — 1 
vanish on the other. 

As we shall see in the next Section, we need the explicit expression for the singlet I'l'^J) and the triplet 
m = 0, ±1, to calculate the ground state of the CBG-Hubbard model. Therefore, we now briefly review 
how to get the projection of the determinantal state |$^^) onto the singlet and the triplet spin subspaces. 

To obtain the singlet |$^^) one has to antisymmetrize each product a||fej|, getting a two-body spin singlet 
operator, and subsequently antisymmetrize with respect to the A'' — 1 indices of the fe^'s; one sees easily that 
this entails the antisymmetrization of the a^'s. Hence 

JV-1 

l<i'^^^)=EH"n<P«|0) (26) 

P !=1 

where a\ ^ creates a two-body singlet state: 

<j-^i<b],-a}^,h\,), (27) 

and P is a permutation of the indices 1, . . . ,N — 1. Since W commutes with the total spin operators and the 
determinant |$^^) is a = state, = as already noted. Applying W on the state in Eq.([26[), one 

can verify this property by direct inspection. 

In a similar way one obtains the triplet projection. We define rf*^ as the two-body triplet creation operator 
with vanishing a-component m = in the orbitals and bj : 

= 7f(4^Ji+«Ii^]T)- (28) 



Then, one has 



io=EH^E-r4n<.«io)- (29) 
p j=i 



^At weak coupling it makes sense to speak about particles in filled shells, which behave much as core electrons 
in atomic physics, and particles in partially filled, or valence, shells. In first order, we must diagonalize the 
perturbation W in the degenerate subspace, dropping the filled shells. Since is a positive semi-definite 
operator, with vanishing eigenvalue, is necessarily the lowest one. Typically, the first-order solution is a 
crude approximation to the exact one, but below, when discussing the exact results, we shall demonstrate that 
many of the qualitative features hold for all U . 
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The m — ±1 components of the triplet in Eq.(|2^) can be otained by means of the total raising and lowering 
spin operators. Let us introduce the following notations 

JV-l 

sf'^ai^.au, s^ = [s+1^ = s-^^f^+^t (30) 

1=1 

JV-l 

st' = bl^b,u s-'^isfW = s-'' = [s+Y (31) 

i = l 

JV-l 

nia=ai^aia, bi = -(rii^ - rii^), n„ = 2_^n,„, n =n|+n|, b = -^{n^ - n^) (32) 



1=1 

JV-l 



it 7 r<zh 1 / " 6 "fa \ "fa X ^ "b "b I '•b rizb "6\ /oo\ 

nia=o'i^bi^, Si = -{riif - Ui^), n^ = ^n,„, n =n^- +n^, S = -{n^ - n^} (33) 

i=l 

(34) 

The states \^]lp^) can be written as 

JV-l 

P i=l 

JV-l 

p i=i 



with 



[: l.m 



Equivalently, the triplet state m = 0, ±1, can also be expressed in terms of the singlet |$^^). It is a 

simple exercise to prove the following identities 



JV-l 



i=i 

JV-l JV-l 



i=i 

JV-l 

E(ftu -"u)l*^p) = E(-"^T +^'T)i*Ap) = (-S^" + = (39) 

i=l 1=1 
JV-l 

J2i-Sr + S-')|-f^° ) = (-S-- + S-")!*^" ) = -V2l$i^-^) . (40) 

i=l 

As for the singlet, one can verify that W\$2p) = 0, m = 0, ±1, using the definitions in Eqs . ( p9| - |3^36[ ) . 

V. THE GROUND STATE AT HALF FILLING 

We are now ready to calculate the ground state |'I'o(Ai', U)) of the half filled CBG-Hubbard model described 
in Section As already observed, gjffjl'i&^j?) ^ good candidate for the non-interacting ground state; its 
quantum numbers are the same as those of |5'o(A'^, U)} and moreover it belongs to the first-order ground state 
multiplet. Using the short-hand notations 

\g°}^g\gl\0), \g''} ^ glgl\0} , [[ggf''^) ^ ±(g]gl - glg\)\0} (41) 

for the three different two-body singlets that one gets from the lowest and the highest energy orbitals g and 
g and 

'•'}^g]g].\0}, |[gg]i-«)^-L(3|g|+pt^t)jo), \[gg]^'-^) ^ glgllO) (42) 
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for the triplet, we propose the following Ansatz for the interacting ground state |\l/o(A'^, U)): 



\<fo{N,U)} = [j,\g''}+T,\g'}] 



l<E>^^)+7o 



E 



(43) 



where the 7's are c-numbers. It is worth to note that in j'I'o(Ai', 17)) the number of particles in the shell Shf is 
a constant given by 2N — 2. A priori, there is no reason for this choice since the total number operator + fi^ 
of particles in the shell Shf does not commute with the Hamiltonian. Nevertheless, we shall see that the 
scatterings which do not preserve this number cancel out provided in Sh f the (2A'' — 2)-body state is a = 
state. We shall refer to this very remarkable property as to the Shell Population Rigidity. We emphasize that 
even constraining the particle- number in Shf to be 2N — 2 with vanishing total spin ^-component, there are 
2N — 2 \ ^ 

^ ^ 1 configurations that can contribute to the ground state expansion in the interacting case, while 

the Ansatz @ contains only |"1>^^) and |<&^"), m = 0,±1. The reason why the W = states |<&^'"-^) with 
S > 1 do not enter into the ground state expansion ( |43| ) comes from the Lieb's theorem the ground state 
must be a singlet and with only two particles outside Shf (in the g and/or g states) the angular momentum 
composition law forbids W = states with 5* > 1. However, we observe that the state [[ffp]"''') <8 I^'aj?) is 
a singlet and |$^^) is a W = state. It has the right quantum numbers and, in principle, it could have a 
non-zero overlap with |^o(A'^, U)}. Nevertheless, the matrix elements of W between this state and the ones in 
the Ansatz are zero. This is the reason why we have dropped it in the expansion (^) . 
The three states in Eq.(p3|) are eigenstates of the kinetic-energy operator Hq: 

Holg") ® |$:i^) = -2t\g") ® 1$^° ) 



1 



(44) 



Expanding c^iCxi occurring in W in a, b, g, g operators of Eq.(^), we may take advantage from the fact 
that \<i>°ll) and \$]ip), m = 0,±1, are W = states. Indeed, since one cannot annihilate g or g over them, 
taking for instance x £ A, 



= c^iC^^^Yf') = E^''-^^'-"'-L"^"tI*af'^) 



(45) 



then, multiplying by pjpj, and the like, one obtains that the contribution proportional to 

^O.ilO'3}\'^ t>{N , U)) 



(46) 



yield nothing. Of course, if a: G B, by the same reasoning we get 

^R^,^RJ,^hlb,^\•^Q{N,U)) = 0. 



(47) 



Hence, we can write 



with 



1 ^ 

^(fli + 3i)(5'T +flT) + 7/= E -^^'^[(ffi +fli)aiT +aa(5T + 9^)] 



(48) 



(49) 



and 



x = N-\-l 



1 ^ ~ ^ 

i^{-9i + gi){-9^ + 9^) + -j^ Y Ri,x-N[{-9i+gi)br-^ + feu(- 



J (50) 
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where :— means that the two sides are equivalent when acting on \'^o{N, U)). Next, we transform the remaining 
two operators, using Eq.(^). We note that the terms containing the sequence g a and similar terms 
with b instead of a and/or some g replaced by g vanish by symmetry; this is evident since Ri,x = for 

any i = 1, . . . , N — 1. Similarly, the sequences like g'' gg also do not contribute. The sequences like a^a^ gg 
are diagonal in the orbital indices of the two operators, since X/^-i Ri,xRj.x ~ Stf, hence, they annihilate 
|^'o(A'^, U)) (recall that in this state all the a orbitals are singly occupied, as we can see by direct inspection 
of the expressions for |$54^) and |"1>^]?), m — 0,±1, of Section [lV[ ). By a particle-hole transformation one 
can show that the remaining terms that create pairs in the shell Shf, like a^a''ga, etc., have no effect, as in 
Eqs.( p^ - |4^ ); we postpone the proof of this last spectacular cancellation until Appendix ^ To sum up, no term 
in W can change the number of particles in Shf, that is, that number is independent of t and 17! This proves 
the remarkable property that we have named Shell Population Rigidity: it holds for any finite A'' and it is not 
specific of the thermodynamic limit. We shall see below that the Shell Population Rigidity characterizes the 
ground state and a suitable subspace of the total Hilbert space. In Fig.^ we have drawn a physical picture 
of the cancellation among diagrams which do not preserve the number Ua + ni,. We can say that the W — 
states entering in our Ansatz are stable with respect to the Hubbard interaction W also in the presence of 
other particles in the system. This fact allows to exactly solve the Schrodinger equation. 

8) is, 8) (a, b) (a, b) (g, g) (a, b) 

+ + + =0 

fe 8) (a, b) (a, b) (a, b) (g, g) (g, g) 

FIG. 2. Cancellation among the scattering amplitudes in the ansatz of Equation ( ^ ) which do not preserve the 
number of particles in the shell Sh, / . 



Taking these cancellations into account one obtains 

U 
AN 



WA~^{g]+9])i9l + gl)igi + gi){9T+g^) + 



and 



+ 2N^ RT(ffi+5l) + (.9T+ffT)"a] +aa(ffT +ffT)] (51) 



Wb := ^(-5t + 9l)i-9l + 9l)i-9l + 3i)(-3T + 9^) + 



AN' 

JV-l 

U 



+ 2ivl^ + [{-9l+9l)b^^+hl{-9^+gr)] ■ (52) 

i = l 

With W written in the transformed picture one can calculate Wl'i/oiN, U)). After very long, but standard 
algebra one finds 

w\g') ® \^r,) = u[^^\9') ® \^r.) + ^\f) ® \^r.)+ 

N-l 

1=1 

JV-l 



2N 

1=1 

JV-l 

2A 



1\ — I 1 

+ W E ^^39]'-'} ® {St" - St'')\<i>TF) \ (53) 



i=l 

and 
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JV-l 

1=1 

JV-l 

+ 2iV E [3!^t(-"^t + '^iT) + 5jffi(jia - "a)] \^°af) + 
1=1 

JV-l ^ 



(54) 



while for the singlet Em=-i(-)"'l[33]^''") ® I^af") one obtains 



N 



E 



(55) 



By means of the identities in Eqs.( p8[ - p9| - }40| ) one can write Eqs.(|53|-p4[) in a compact form 



27V 



and 



Wig") ® I*"-" ) = U ■ 



2iV-l,_o 



2iV 



(56) 



(57) 



Let us now consider the second row in Eq.(^5|). If our Ansatz is correct, the state in the square brackets 
must be proportional to |$^°). In this way one could close the equations and find an exact eigenstate. In 
Appendix |^ we prove that this is the case and in particular that 

JV-l 



E 

1=1 

Hence, Eq.(|55|) yields 

1 



V2{N' - 1)|<1>^° ) 



(58) 



jV + 1 
iV 



(59) 



This result, together with Eqs.( |5^ - |5^ ) and Eqs.(^), allows us to reduce the Schrodinger equation (i^Hubbard- 
E)\'io{N, U)} = to the diagonalization of the matrix 



H{N,U) 



-I- 2JV 2JV 



u 

2JV 



-t- 2JV 2JV ^ 



V2(JV^-1) , V2(JV^-1) , 

2JV ^ 2JV ^ 



(60) 



where the relations 



m ® (K^A^i] Op") ® = m ® [\f) ® 11-:^^)] 



(61) 



and 
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m— — 1 m— — 1 

1 1 

m — — 1 m— ~1 

have been used to express X/m=-i ("^af ) terms of {^'afI^'af) ■ 

The thermodynamic limit A'^ — > oo is weU defined and non-trivial if we rescale t in such a way that Nt = 
t =const (which corresponds to the case of a fixed energy gap between the first two energy levels). The 
Hamiltonian matrix in Eq.(^) becomes 

/ -2t + U U/V2 
H(oo,U)= 2t + U ~U/V2 I (63) 

y U/V2 U 

with eigenvalues 



So^U, £± = U±A, A= + 4*2 . (64) 
The (unnormalized) eigenvector corresponding to the lowest eigenvalue f _ gives 



1 



|*o(cx),C/)) = 



(t - + ^)) Is") + ^1/)] ® |4>:i» ) + V2(2 + ^) E (-)™|[55]^"") ® l<i>i^.'") . 

m=-l (65) 

So far we have proved that the Ansatz in Eq.(^) gives three exact eigenstates of the CBG- Hubbard model 
at half filling. The eigenvalue of lowest energy £- reduces to the non-interacting ground state energy for 
f/ — > 0, while the associated eigenstate (^) to the state Ig") ® |$^^). On the other hand, when U oo the 
original Hamiltonian can be mapped onto what we can define the CBG-Heisenberg Hamiltonian: 

-ffuoiscnborg = — ^ E(^!= ' ^j; - ^) . (66) 

xeA yes 

This model is exactly solvable and the ground state is obtained by projecting onto the singlet the state where 
particles lie on the ^-sites with spin up and the remaining A'' on the Z3-sites with spin down (Neel state): 
one finds 

N 

|*o(iV, ^)) = n -I.P(.+iv)!0) (67) 



p 



with al y — -^{^ti'^li ~ '^I'l'^vt'^' ground state energy is —2t^/U for N ^ oo and is equal to the first order 
approximation of £- in the small parameter (t/U)^ (the same conclusion holds for any finite A'^, although it 
is more tedious to prove). Furthermore, by direct inspection of the lowest energy eigenvector of H{N, U) one 
can show that |^o(-'V, U)) becomes the state in Eq.(|67[) when U ~-> oo. Therefore, |^o(-^5 U)) reduces to the 
ground state of the CBG-Hubbard model in the two opposite limits U ^ and U ^ oo. 

To prove that the lowest energy eigenstate of H{N, U) is the unique ground state we are looking for, one can 
exploit the ground state uniqueness of the Heisenberg Hamiltonian proved in Ref. Since |\['o(A'', oo)) is the 
ground state of the Heisenberg model and since the ground state of the half-filled Hubbard model is unique, 
a level crossing for some value of U would be required if |'I'o(A', 17)) were an excited state, contradicting the 
uniqueness. 

In conclusion we have proved that the Ansatz in Eq. (^) with (7g ,79,70) given by the lowest energy eigen- 
vector of the matrix in Eq.(^^ is the exact ground state of the half-filled Hubbard model defined on the CBG 
of size |A| = 2A' and repulsion U . In the next Section we shall calculate some physical quantities, as the 
spin-spin correlation function, and we shall prove that the particles are antiferromagnetically ordered. 



VI. RESULTS AND DISCUSSION 



The exact ground state solution for arbitrary but finite TV allows to study the ground state energy Eo{N, U) 
as a function of A' and U. Taking tN = t = 1, in Fig.(^) we have plotted Eo{N, U) in the range A' = 1, .., 10 
and < t/ < 20. 
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E,(N,U) 




FIG. 3. Ground state energy Eq{N, U) of the half-filled Hubbard model on the CBG as a function of the on-site 
repulsion parameter U and the number of lattice sites A'^ in each sublattice. The hopping parameter has been chosen 
to be t = 1. 



The figure shows that Eo{N, U) is a monotonically increasing function of A'^ and U, that is Eo{N + 1, U) 
Eo(N) > and dEo{N, U)/dU > 0. The limit iV oo at fixed U is given by £^ in Eq.(||): 



lim lim E'-'^\\A\,t,U) = U- ^/U^+Af^ < (68) 

where i?^'^'(|A], t, U) is the internal energy of the system at inverse Temperature /3 and number of lattice sites 
|A| = 2N. On the other hand, one can calculate the internal energy at zero Temperature in thermodynamic 
limit; E'-'^\\A\ oo, t,U) = + ksS/fJ + ^|A|n with ^0 from Eq.(|l|) and 5 the entropy. Following the 
procedure of van Dongen and VoUhardt one gets 

!7{;3)/i Ai w^^ ( -2t = -2Ni < if Hs not rescaled , , 

iT^IA^^^oc^ (1^1'*'^) = I 0(1) iff- 1/|A| ■ (69) 

In the case t=const the two limits commute and the trivial thermodynamic behavior of the model originates 
from the infinite gap between the lowest level and the zero-energy levels: the two particles in the g orbitals are 
completely decoupled from the dynamics of the system. On the other hand, if t is rescaled the gap remains 
frozen and taking |A| — > oo first, we can only say that the ground state energy is 0(1), but we cannot predict 
the exact amount within the scheme proposed by van Dongen and VoUhardt. This is a consequence of the 
fact that in the limit |A| —> oo, we retain only the extensive contributions to the internal energy. 

The ground state energy Eq.(p^) is always higher than that of the non- interacting case and hence non-trivial 
correlations survive in the thermodynamic limit. We have calculated the ground state average of the number 
of doubly occupied sites: 

— W(N U) 1 1 

D{N, U) = ' = ■^{*o(iV, t/)! W|*o(iV, 17)) = ■^{*o(A^, ?7)|//Hubbard - i?o|*o(iV, U)) = 



li5o(iV, U) + |[7s(^, U) - jUn, U)] (70) 



where 73 and jg are the first two components of the normalized ground state vector, see Eq.(p^). In Fig.(|ja) 
we have plotted the trend of D{N,U) as A*" increases for different values of U. As expected D{N,U) is a 
monotonically decreasing function of U, dD{N,U)/dU < 0, approaching zero for [/ — > 00 [where the exact 
ground state reduces to the antiferromagnetic Neel state, see Eq.(^^]. Nevertheless, D{N,U) shows a non- 
trivial behaviour as A'" increases at fixed U values. In the weak coupling regime, U <^t, the number of doubly 
occupied sites grows as A'^ becomes larger and larger converging to a finite value when A' — > 00. The opposite 
trend is observed in the strong coupling regime, U ^ t, where D{N, U) decreases as A'' increases. In the 
intermediate regime U t, D(N, U) is an increasing function of A'^ for small A', but becomes a decreasing 
function for large A''. Hence, there is a critical value Uc{N) where the analytic continuation of D{N,U) to 
real A'^ verifies 



dUdN J 



= (71) 

(iV) 
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D W 




FIG. 4. a) Trend of the ground state average of the number of doubly occupied sites versus N in the range 
1 < Af < 20 for different values of the Hubbard interaction parameter U =0.1, 1, 2, 3, 5, 10. h) Ground state average of 
the interaction term W as sl function of U in the range < !7 < 10 for three different number of sublattice sites, A'^ = 
1, 2, 10. The hopping parameter has been chosen to be i = 1 in both cases. 



In Fig.(y6) we report the ground state average of the interacting term W as a function of U for three 
different values of A'^. There is always a critical value for U where the repulsion is maximized. Moreover, 
limjv^oo W{N, U) = Woo{U) 7^ for any finite U and the system cannot avoid double occupation neither in 
the thermodynamic limit. In the limit U — » 00, W{N, U) ^ jj, that is, W{N, [/) — > in agreement with the 
Gutzwiller Ansatz psf . 

Next, we have calculated the spin-spin correlation function 

Gspir.[x,y) = {<^o{N,U)\SlSl\^o{N,U)) = lim Gspi„(a:,2/;/3,/i = 0, |A|) = lim lim Gspin(a;,2/; /i, 1A|) 

/3^oo /3— ^00 /i^O (72) 

where Gspin{x, y; I3,h, |A|) was defined in Section S| — ^{uxt; — fixi) and |\E'o) is the normalized ground 
state 



1*0) 



m — — 1 



(73) 



with AfS = {^afI^af), a/"" = (^afI^af), rn = 0,±1, and 79 + 7s + 7o = 1- Due to the Sjv (g) Sjv ® Z2 
symmetry, Gspin{x,y) can be written as 



{Go X = y 
Gon X e A{x e B) a.ndy e A {y e B) 
Goff X e A {x e B) and y e B [x e A) 



(74) 



and by exploiting the sum rule '^y^f^GBpin{x,y) — Go + (N — l)Gon + NGoB ~ 0, one can express Gon in 
terms of Goff and Go: 



Go 



N 



-Go 



Go 



TV - 1 N-1 
The problem is then reduced to the evaluation of Go and Goff. We have 



(75) 



Go = G.pin{x,x) = {*oi(5|)2i*o) = ^ ^(vPoKSD^lvI/o) = 



(76) 



More effort is needed to calculate Goff. Expanding the c operators in g, g, a and b operators of Eq.(|7|) and 
taking into account that (X^^-g^ 5'|)^j^'o) = 0, we get 

^°« = -i^EE<*oi^^^^i*°> = -i^ j^(*oi ('^(-)'^(ff:+ff:)(3.+3.)') 1*0) + 

+ 5(*olS^''^(-)^(gI +5l)(5. +ff.)l*o) + (^olS^-S^^l^o) (77) 
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with {—)" — for a =|, j. After long but standard algebra one obtains 



GoS ~ 



;^<g [(79-79) +270J-Y3 ^^^^^^^^^ {'^afI^ I^Af) + ("^AfP i'l'AF) + 



+ 



{$Ji-V^11'.^F'>-(*iFl5^1*iF)+ E (^AFl^^^S^^l^i^F) 



and exploiting the Z2 symmetry —^ hi and 6i — > a^, which implies — > (— )'^1'I'af*^ )> 



(^afIs^I^af) = i(*iFl^^" + s^''l*iF) = 



(78) 

(79) 
(80) 



(81) 



where in Eg. (|79[) we have used the identity Eq.(|39|). Finally we recall that the W — Q states |$^'^^^) are 
eigenstates of the square of the total spin operators of each sublattice S\ = (S*^")^ + \ {S'^°'S~°' + S~°'S^°') 
and Si = {S'^'f + i(S+''5'-' + S-''S+'') with eigenvalue (i^)(i^ + 1) = \{N^ ~ 1); hence 



Ms = -S 

Substitution of Eqs.(|7|-|ol-|l[||) in Eq.(||) yields 



A' -1 Ir, n2 2, 1 2 /iV2_l 

^2 + 8 [(79 - 79) + 270I - 370 - Y -^g 70(79 - 79) 



(82) 



(83) 



where we have taken into account Eq. ( p2[ ) to evaluate the ratio Mi /No = \J {N^ — 1) /3. 



a) 



1 G'off 



ft) 



N^lO 

N=3 

N=2 



-0.075 
-0.1 
-0.125 
-0 . 15 



u 



N=10 
N=2 



FIG. 5. a) Gon versus U in the range < f/ < 20 for three different vaiues of the number of sites N — 2, 3, 10. 
b) Goff versus U in the range < f7 < 20 for three different values of the number of sites N = 2, 3, 10. The hopping 
parameter has been chosen to be t = 1 in both cases. 



In Fig.(|5|) we report the trend of Gon and GoB versus U for three different values of = 2, 3, 10. According 
to the Shen-Qiu-Tian theorem Gon is always larger than zero while Goff is always negative. Now consider 
the ground state average of the square of the staggered magnetization operator 

m\p = j^{^o\lY,<^)Slf\'^o), e{x) = l, -lior x€A, B . (84) 
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The Shen-Qiu-Tian theorem imphes that each term in the expansion of Eq.(p4|) is non-negative. We emphasize, 
however, that for |^| — \B\ this does not imply that m\p is an extensive quantity! Remarkably, our results 
show that m\jr — Go + {N — l)G'on — NGoB is extensive for any value of the on-site repulsion parameter U 
and constitute the first example of antiferromagnetism in the ground state of an itinerant-electrons model of 
arbitrary size. 

We also observe that the thermodynamic limit yields a non-vanishing result: 

lim Go = 7, lim Gon = - lim GoS = (85) 

|A|-.oo 4 |A| 

or 

lim lim lim |Gspin(2;, y; /3, /i, |A|)j = for a; / y . (86) 

A|— >oo /3^oo /i^O iZ 

By comparing Eq.(p6|) with Eq.(p^), we see that the lim/3_,oo limh^o do not commute with the lim|A|^tx>- 



VII. SUMMARY AND CONCLUSIONS 



We have introduced the CBG-Hubbard model and explicitly written the ground state wave function at half 
filling for any repulsion parameter U . The topology of the graph A allows (27V — 2)-body eigenstates \'^^j^p'^) 
of the Hamiltonian which are free of double occupation. They can be obtained projecting the determinantal 
state |$a_f) on a given spin-S subspace. |$af) has the antiferromagnetic property, that is the map A ^ B is 
equivalent to a spin-flip. 

We have shown that the ground state has a fixed number of particles in the zero-energy one-body shell Shf, 
independent of t and U. This very remarkable Shell Population Rigidity holds in any finite-size system and 
it provides the closure of the equations. This implies extra-conservation laws in a suitable subspace of the 
full Hilbert space including the ground state. Only the singlet l"!?^^) and the triplet ), m = 0,±1, are 

involved in the ground state expansion. Qualitatively, we may say that according to Eq.(g5|) the particles in 
the shell Shf manage to avoid the double occupation. 

We have calculated the spin-spin correlation function and shown that the ground state exhibits an antifer- 
romagnetic order for any non-zero U even in the thermodynamic limit. The kinetic term induces non-trivial 
correlations among the particles and an antiparallel spin configuration in the two sublattices turns out to 
be energetically favoured. Therefore, the model contains the basic ingredients to understand how the subtle 
competition between the delocalization induced by Ho and the localization induced by W may give rise to a 
magnetically ordered state even outside the strong coupling regime (where the Hubbard model is equivalent 
to the Heisenberg model). On the other hand, within the scheme proposed by van Dongen and VoUhardt, 
the model is well described by a paramagnetic Hamiltonian; the difference stems from taking the thermody- 
namic limit before the limit of zero Temperature, because then the kinetic term completly decouples from the 
interaction and loses any role in the physical response functions. 

The present formalism lends itself to solve more realistic Hubbard Hamiltonians with an increasing number 
of negative and positive energy levels. Higher-spin projections \^^p^} are involved in the ground state 
expansions; the results will be published elsewhere. We are currently investigating group-theory aspects of 
this model and its extensions. 

Finally, we recall [p3|[ that the standard Hubbard model on a A'' x square lattice and periodic 

boundary conditions also has |$^'^^^)-like eigenstates. It could be that the Shell Population Rigidity holds in 
this case too, but the proof is lacking. 



APPENDIX A: SHELL POPULATION RIGIDITY 

We have seen in Section ^ that the most part of the contributions which do not preserve the number of 
particles in the shell Sh / yield nothing. To prove the Shell Popoulation Rigidity property we need to show that 
the sequences like a^a^ga and similar terms with 6 instead of a and/or g replaced by g vanish or in formula 

R,,.Rj,'.Rm,^al^al^ama\^AF'') = (Al) 

ijm 
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and the like with a replaced by b. This last spectacular cancellation can be seen by particle- hole transforming 
Eq.(Al). Under a particle-hole transformation Cx<t — > cl^a and the W — Q states |$^'^^^) Is'^s'') ® I'^'af*^^)) 
modulo a phase factor. Hence, Eq.(|Al|) is equivalent to 



x^A. ijm 



(A2) 



Let the spin of the operator be up for the sake of definiteness (the same reasoning holds in the case of 
down spin). The l.h.s. of Eq.(A2) can be rewritten as 



ijm j 



(A3) 



where we have exploited the lack of double occupation of the W = state, see Eq.(^). Therefore, we are left 
to prove that 



(A4) 



We recall that R is an (A'^— 1) x rectangular matrix whose rows are orthonormal vectors which are orthogonal 
to the A^— dimensional vector (1, 1, . . . , 1, 1). It is a simple exercise to verify that 



R 



1 


1 





V2 


V2 


1 


1 


2 


V6 






1 


1 


1 



\ ^N(N+l) ^N(N+-L) ^N{N+l) i/]V(JV+l) 



(iV+1) / 



x> m + l 

Rm,x = { x = m + l 

-a/^ a;<m (A5) 



is a correct choice and that Eq.(A4) is identically verified. 



APPENDIX B: PROOF OF EQ.(M) 



To prove Eq.(|58|) we shall use the definitions (gg[j35|-g6|) for the triplet W = state \^2p), m = 0, ±1. The 
first term yields 



fc p k j 

= E(-)"E(^"-^J[;,)^^Sn<p«io) + 

+E(-)"E^^SE(^r-^J«vU) n 4..) 

= ^E(-)"E-U)n<.«io) + 

-^^EH"EE^^S^:^^)^<..)!o)- 



P 3 kytj i^],k 

V2(;v-i)|$-)-y2E(-)^EE^.^iW^^^) n <pi. 

P j ki^j i7^j,k 



(Bl) 



while second term yields 



17 



^(n^l -rifcj, - "-fcT +"-fci)l*AF) = --r'ti-np(k)r +fh'(k)i)'^ 

k P k j i^j 



+ ^(-)''^'^tp(3)^^'^"k^ -^H-np^k)t +np^k)iWlp(^k) n '^Ip(>)|0) 

P J k^J 

n <PH^o^ 



+ 2E(-)^EE 



j,P{j)'k.P(k) 11 ^i,P(i)l 

(0)t ^(0)t 



(B2) 



P i ki^i ii'j,k 

The third term can be computed following the same steps of Eq.(p3l|) and the final result is 

^(+i)t (_i)t 
j,p(j)'fe.p(fe) 

* P 3 k^j ij^j,k 



E(-sr + = V2iN - mrp) - ^E(-)" EE 11 <Pi^^\^ 



(B3) 



Denoting by [..] the left hand side of Eq.(|58[), from Eqs.(Bl-B2-B3) one obtains 



i^j.k 



(B4) 



Let us consider the first and the third term in the square bracket of the right hand side of Eq.(B4). We have 

Z^'- [j,p'U) k,p'{k)+^j,p'uyk,p'{k)\ 11 '^»,p'{»)iu/ - 

P' iT^i.fc 

~ > yj,p(k)^k.p(j) + ^],p{k)^k,p{j)\ 11 '^i,p(!)i'^'' 



ii'3,k 

where the permutations P and P' satisfy 

P'{k)^P{j), P'{j) = P{k), P'{i) = P{i) yij^j,k; 

^P' 



hence (- 



(— ) . Substituting Eq.(B5) into Eq.(B4) and taking into account that 



j-ry'(+iy , ^(o)t ^(o)t 



'j,p{k)'k,p{j) ^ 'j,p(j)'k,p{k) ^ ' j,p{k)'k,PU) - j,PU)'^k,p{k) 
one obtains 

[..] = 3V2(iV- 1)1$^^°) + V2E(-)'^EE4pU)'^I,P(fc) n = ^(^' - l)l*°ip) 



i¥'j,k 



(B5) 



(B6) 



(B7) 



(B8) 



that is, Eq.i 
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